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ABSTRACT 


The brachyuran crab Bythograeidae Williams, 1980 is common in hydrothermal vent fields worldwide and has 


recorded to sixteen species of six genera. In this study, we firstly determined the cytochrome c oxidase subunit 1 
(CO1) DNA barcodes for the fifth species of Austinograea, A. hourdezi, from hydrothermal vent regions of the North 
Fiji Basin in southwestern Pacific Ocean. All CO7 DNA barcodes of A. hourdezi were identical. The interspecies 
variations of three bythograeid genera were 10.9-13.3% for Austinograea, 6.6-15.7% for Bythograea, and 9.7% for 
Gandalfus. These results would be helpful to understand taxonomy of brachyuran crabs living in hydrothermal vent 


fields using CO7 DNA barcodes. 
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INTRODUCTION 


The brachyuran crab Bythograeidae Williams, 1980 is com- 
mon in hydrothermal vent fields worldwide and has record- 
ed to sixteen species of six genera (Guinot and Segonzac, 
2018). For DNA-based species identification, it is generally 
considered a suitable DNA barcode when sequence varia- 
tion among species is 10 times greater than those within 
species (Hebert et al., 2003; Matzen da Silva et al., 2011). 
Based on these criteria, the mitochondrial cytochrome c ox- 
idase subunit 1 (CO7) DNA barcodes for bythograeid crabs 
were confirmed worth supporting the traditional taxonomy 
and utility recognizing the existence of sibling (or cryptic) 
species (Mateos et al., 2012; Kim et al., 2013; Guinot and 
Segonzac, 2018). 

Austinograea hourdezi Guinot & Segonzac, 2018 is found 
across wide hydrothermal vent areas in the southwestern Pa- 
cific Ocean (Guinot, 1989; Guinot and Segonzac, 2018) (Fig. 
1A). It is distinguished from its relatives with spots on outer 
surface of chelae in males, patches of setae on inner surface 


of chelae, subhepatic region and meri, and ratio between the 
Ist and 2nd pair of gonopods (Guinot and Segonzac, 2018) 
(Fig. 1B-D). In this study, we determined firstly the CO7 
DNA barcodes of A. hourdezi collected from the hydrother- 
mal vents of the North Fiji Basin in southwestern Pacific 
Ocean. Additionally, we obtained mitochondrial 16S rDNA 
sequences to compare the 16S rDNA sequence of A. hour- 
dezi retrieved from GenBank. 


RESULTS AND DISCUSSION 


In November 2016, six A. hourdezi individuals were collect- 
ed from three hydrothermal vent regions (1 ind. from 
18°49'S, 173°30'E, 2,719 m depth; 3 inds. from 18°51'S, 
173°30’E, 2,721 m depth; 2 inds. from 19°03'S, 173°29'E, 
2,717 m depth) (Fig. 1A) in the North Fiji Basin. Morpho- 
logical identification, genomic DNA extraction, sequencing, 
sequence annotation, and nucleotide diversity analyses fol- 
lowed the methods of Guinot and Segonzac (2018) and Kim 
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CO1 DNA Barcoding for Bythograeid Crabs 


Fig. 1. Photograph of the hydrothermal vent chimney in the North Fiji Basin of the Southwestern Pacific Ocean where Austinograea 
hourdezi Guinot & Segonzac, 2018 was collected. A, Overview of sampling site and the surrounding community of species including 
bythograeid crabs Austinograea spp., alvinocaridid shrimp species, provannid snails A/viniconcha boucheti Johnson et al., 2015 and 
A. kojimai Johnson et al., 2015, mytilid mussel Bathymodiolus sp., chiridotid sea cucumber Chiridota sp.; B-D, Diagnostic keys of A. 
hourdezi based on Guinot and Segonzac (2018). Arrowheads indicate patches of setae on inner surface of chelae (B), spots on outer 
surface of chelae in males (C), and patches of setae on subhepatic region and meri (D). Scale bar=20 cm. 


Table 1. Sequence divergence among the partial sequences of mitochondrial 16S rDNA and cytochrome c oxidase subunit 1 (CO1) 
genes from Bythograeidae. Nucleotide sequence variations (%) were calculated using the p-distance method in Mega 10.0.5. Upper 
right is for 16S and lower left for CO1. The bold letters indicate interspecies variation within the genera Austinograea, Bythograea and 
Gandalfus 


1 

(7°, 0.22%) 2 3 4 5 6 7 8 9 10 1i 12 13 14 
1. Austinograea hourdezi (6°, 0.00”) = 7.4 6.1 6.3 10.4 14.9 15.1 14.4 14.5 15.3 93 98 9.9 13.2 
2. Austinograea alayseae 10.9 = 5.6 5.8 12.2 15.6 14.4 15.1 14.1 15.1 10.2 90 9.4 13.6 
3. Austinograea rodriguezensis 11.7 11.1 - 4.5 10.0 13.1 13.1 13.6 13.0 13.8 7.3 7.3 7.4 11.6 
4. Austinograea williamsi 11.1 10.3 13.3 - 10.0 13.3 13.3 13.0 13.2 14.0 86 86 8.1 12.9 
5. Allograea tomentosa 15.9 15.7 13.3 15.9 - 13.0 13.0 13.5 13.7 13.2 7.3 9.6 8.8 11.6 
6. Bythograea galapagensis 17.7 15:7 16:5 17.1 15.1 - 7.3 5.7 4.7 7.1 11.3 12.2 11.9 12.4 
7. Bythograea laubieri 16.3 15.3 15.5 16.1 14.9 13.9 - 7.1 8.8 0.7 11.7 11.7 12.1 13.0 
8. Bythograea microps 17.5 18.5 18.3 18.7 17.9 15.5 15.7 - 8.6 7.3 11.9 13.8 12.6 13.5 
9. Bythograea thermydron 19.3 17.3 19.7 17.9 18.3 9.3 14.1 15.1 - 9.2 12.1 11.3 11.3 12.5 
10. Bythograea vrijenhoeki 16.1 15.7 16.5 16.3 15.7 15.1 6.6 15.1 14.5 - 11.9 12.0 12.3 13.2 
11. Cyanagraea praedator 17.3 17.3 15.5 18.3 12.3 14.1 15.9 18.1 15.5 16.7 - 96 88 9.1 
12. Gandalfus puia 13.1 14.1 14.7 14.9 12.7 16.7 14.1 16.3 16.7 13.9 14.5 - 5.2 9.8 
13. Gandalfus yunohana 15.9 14.9 15.7 14.3 13.5 16.7 16.1 18.3 17.9 16.3 15.9 97 - 9.2 
14. Segonzacia mesatlantica 16.5 15.9 15.1 15.9 14.3 15.5 14.9 15.7 18.5 15.5 13.5 14.3 14.9 - 


The accession No. of bythograeid crabs are A. hourdezi, J1Q407451, MK382751-MK382753, MK359095-MK359097, MK359098-MK359103; A. alayseae, 
KC851803 (mitogenome); A. rodriguezensis, 1Q407452 and JQ407476; A. williamsi, J]Q407448 and JQ407474; A. tomentosa, JQ407446 and JQ407472; 
B. galapagensis, 1Q407444 and JQ407470; B. laubieri, 1Q407442 and JQ407468; B. microps, JQ407441 and JQ407467; B. thermydron, JQ407440 and 
JQ407466; B. vrijenhoeki, 1Q407443 and JQ407469; C. praedator, JQ407445 and JQ407471; G. puia, NC_027414 (mitogenome); G. yunohana, NC_013713 
(mitogenome); S. mesatlantica, 1Q407447 and JQ407473. 

*Number of individuals. 

>Intraspecies variation (%). 
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et al.(2013). All specimens and genomic DNA extracts of A. 
hourdezi are deposited in Korea Research Institute of Bio- 
science and Biotechnology. 

The newly obtained sequences for A. hourdezi were regis- 
tered to GenBank (accession No. MK382751-MK382753, 
MK359095-MK359097 for 16S rDNA; MK359098-MK 
359103 for CO/) (Table 1). The lengths of alignments were 
446 bp for 16S rDNA and 503 bp for CO1. All CO1 sequenc- 
es were identical, while 16S rDNA sequences detected a 
variable nucleotide site in the single individual (MK382752). 
Interspecies variations of 16S rDNA and CO/ within the ge- 
nus Austinograea were 4.5-7.4% and 10.3-13.3%, respec- 
tively. In the case of the genus Bythograea, interspecies vari- 
ations were 0.7-9.2% for 16S rDNA and 6.6-15.7% for CO1. 
Two Gandalfus species showed the variations of 5.2% for 
16S rDNA and 9.7% for CO1. The sequence diversity of CO7 
was two times more than that of 16S rDNA. Meanwhile, the 
sequence variation of 16S rDNA between B. laubieri and B. 
vrijenhoeki was only 0.7%, which is considered to represent 
the most recently derived lineages or cryptic species. 

Recently, A. aff. williamsi collected in the southwestern 
Pacific Ocean was revised as the fifth species of Austinograea, 
A. hourdezi, by Guinot and Segonzac (2018). As of 8 Jan. 
2019, GenBank contains the sequences of mitochondrial 16S 
rDNA and cytochrome b, and nuclear 28S rDNA from one 
single individual of A. aff. williamsi (JQ407451, JQ407421, 
and JQ407435), but no sequence for A. hourdezi. These re- 
sults would be helpful to understand taxonomy of brachyu- 
ran crabs living in hydrothermal vent fields using COJ DNA 
barcodes. 
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